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ABSTRACT 

We identify 4 unusually bright (Hi6o,ab< 25.5) galaxies from HST and Spitzer CANDELS data with 
probable redshifts z^7-9. These identifications include the brightest-known galaxies to date at z > 7.5. 

As F-band observations are not available over the full CANDELS program to perform a standard 
Lyman-break selection of z > 7 galaxies, we employ an alternate strategy using deep Spitzer/IRAC 
data. We identify z~7.1-9.I galaxies by selecting z>6 galaxies from the HST CANDELS data that 
show quite red IRAC [3.6] —[4.5] colors, indicating strong [OHI]-fH/3 lines in the 4.5/im band. This 
selection strategy was validated using a modest sample for which we have deep Y-band coverage, and 
subsequently used to select the brightest z >7 sources. Applying the IRAC criteria to all HST-selected 
optical-dropout galaxies over the full ~900 arcmin^ of the CANDELS survey revealed four unusually 
bright z ^ 7.1, 7.6, 7.9 and 8.6 candidates. The median [3.6] — [4.5] color of our selected z ~ 7.1-9.1 
sample is consistent with rest-frame [OHI]-|-H/3 EWs of ~1500A in the [4.5] band. Keck/MOSFIRE 
spectroscopy has been independently reported for two of our selected sources, showing Lya at redshifts 
of 7.7302±0.0006 and 8.683to.oo4, respectively. We present similar Keck/MOSFIRE spectroscopy for 
a third selected galaxy with a probable 4.7cr Lya line at Zgpec =7.4770±0.0008. All three have Higo- 
band magnitudes of ~25 mag and are ~0.5 mag more luminous (Migoo ~ —22.0) than any previously 
discovered z^8 galaxy, with important implications for the UV LF. Our 3 brightest, highest redshift 
z>7 galaxies all lie within the CANDELS ECS field, providing a dramatic illustration of the potential 
impact of field-to-field variance. 

Subject headings: galaxies: evolution — galaxies: high-redshift 


1. INTRODUCTION 

The first galaxies are believed to have formed within 
the first 300-400 Myr of the Universe and great strides 
have been made towards identifying objects within 
this era. Since the installation of the Wide Field 
Camera 3 (WFC3) instrument on the Hubble Space 
Telescope (HST), an increasing number of candidates 
have been identified by means of their photometric 
propertie s, with >700 probab le galaxies identified at 
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One of the most interesting questions to investigate 
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with these large samples is the build-up and evolution of 
galaxies. While these issues have long been explored in 
the context of fainter galaxies through the evolution of 
the UV LF, less progress has been made in the study of 
the most luminous galaxies due to the large volumes that 
must be probed to effectively quantify their evolution. 

The entire enterprise of finding especially bright galax¬ 
ies at z > 7 has been limited by the availability of suf¬ 
ficiently deep, multi-wavelength near-infrared data over 
wide areas of the sky. The most noteworthy such data 
sets are the UKIDSS UDS program (Lawrence et al. 
2007), the UltraVISTA program (McCracken et al. 2012), 
the 902-orbit CANDELS program from the Hubble Space 
Telescope (Grogin et al. 2011; Koekemoer et al. 2011), 
the BoRG/HIPPIES pure-parallel data set (Trent! et al. 
2011; Yan et al. 2011; Bradley et al. 2012; Schmidt et al. 
2014; Trent! 2014), and the ZFOURGE data set (Tilvi 
et al. 2013; 1. Labbe et al. 2016, in prep) 

Of these surveys, arguably the program with the best 
prospects for probing the bright end of the z > 7 pop¬ 
ulation would be the wide-area CANDELS program|3 
The challenge with CANDELS has been that it is only 
covered with particularly deep near-infrared observa¬ 
tions from 1.2^m to 1.6^m, but lacks HST-depth Y- 
band observations at l.OS^m over the majority of the 


® In principle, the wide-area (~1 deg^) UDS and UltraVISTA 
programs have great potential to find large numbers of bright z'>6> 
sources as demonstrated by the recent Bowler et al. (2014) results 
(see also Bowler et al. 2015), but may not yet probe deep enough 
to sample the 2 : > 8 galaxy population. 
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area. Deep observations at 1.05/rm are needed for 
the determination of photometric redshifts for galax¬ 
ies in the redshift range z ~ 6.3 to z ~ 8.5. While 
this can be partially compensated for by the availabil¬ 
ity of moderate-depth 1.05/rm observations from var¬ 
ious ground-based programs over the C ANDELS pro- 
gram, e.g., HUGS (iFontana et al.l 1201411 , UltraVISTA 
(|McCracken et al.l 1201211 . and ZFOURGE (1. Labbe et 
al. 2016, in prep), such observations are not available 
over the entire program, making it difficult to consider a 
search for bright z > 7 galaxies over the full area. 

Fortunately, there appears to be one attractive, alter¬ 
nate means for making use of the full CANDELS area to 
search for bright z > 7 galaxies. This is to exploit the 
availability of uniformly deep Spitzer/IRAC observations 
over the full area (e.g., Ashby et al. 2013) and redshift 
information present in the [3.6] —[4.5] colors of 2 ~ 5- 
8 galaxies. As demonstrated by many authors (e.g., 
Labbe et al.l 120131: iSmit et al.l 12014 l2015t [Bowler et a.lJ 

20l4lLar)orte et al.ll2014aLl2015HHuang et al.ll2016ll . the 

[3.6] —[4.5] colors appear to depend on redshift in a par¬ 
ticularly well-defined way, a dependence which appears 
to arise from very strong nebular emission lines, such as 
Ha and [0111]A5007 A, which pass through the IRAC 
bands at particular redshifts. For example, while z ^6.8 
galaxies have very blue [3.6] —[4.5] colors likely due to 
contamination of the [3.6] filter by [0111]-|-H/3 lines (and 
no similar contamination of the [4.5] band), z >7 galaxies 
exhibit much redder [3.6] —[4.5] colors, as only the 4.5^m 
band is contaminated by the especially strong [0111]-l-H/3 
lines ([Labbe et al.ll2013i : IWilkins et al.ll2014 ISmit et ahl 

IWl . 

Here, we make use of the redshift information in the 
Spitzer/IRAC observations and apply a consistent set of 
selection criteria to search for bright z ^ 8 galaxies over 
all 5 CANDELS fields. A full analysis of the HST -1- 
ground-based observations is made in preselecting can¬ 
didate z > 6 galaxies, for further consideration with the 
available Spitzer/IRAC data. The identification of such 
bright sources allows us to better map out the bright end 
of the UV luminosity function (LE) at z > 7 and con¬ 
strain quantities like the characteristic lu minosity M* or 
the fu nctional form of the LE at z > 7. iBouwens et al.1 
(|2015fl only observe a modest (^0.6±0.3 mag) brighten¬ 
ing in the characteristic luminosity M* - or bright end 
cut-off - from z ^ 8 to z ~ 5 taking advantage of the full 
CANDELS -b XDF -b HUDF09-Ps search area (-1000 
arcmin^). ? also report evidence for a limited evolution 
in the characteristic luminosity with cosmic time, based 
on a wider-area search for z — 6-7 galaxies found over 
the —1.7 deg^ UltraVISTA-bUDS area. Limited evolu¬ 
tion was also reported by Finkelstein et al. (2015) in 
subsequent work, but utilizing a —3-15x smaller area 
than Bouwens et al. (2015) or Bowler et al. (2015) had 
used. 

This paper is organised as follows: §2 presents our 
z —5-8 catalogs and data sets, as well as methodol¬ 
ogy for performing photometry. §3 describes the selec¬ 
tion criteria we define for our samples and methodol¬ 
ogy. §4 presents the results of our investigation and dis¬ 
cusses the constraints added by Y -band observations and 
Keck/MOSEIRE spectroscopy. In §5, we use the present 
search results to set a constraint on the bright end of 


TABLE 1 

Summary of Data Sets Utilized in Current Search. 


Data Set 

Area 

J 125 

Depth (Sc) 
Hiqo [3.6] 

[4.5] 

CANDELS GS DEEP 

64.5 

27.8 

27.5 

26.1 

25.9 

GANDELS GS WIDE 

34.2 

27.1 

26.8 

26.1 

25.9 

ERS 

40.5 

27.6 

27.4 

26.1 

25.9 

GS other® 

31.8 





CANDELS GN DEEP 

62.9 

27.7 

27.5 

26.1 

25.9 

GANDELS GN WIDE 

60.9 

26.8 

26.7 

26.1 

25.9 

GN other® 

34.0 





CANDELS UDS 

191.2’" 

26.6 

26.8 

25.5 

25.3 

CANDELS COSMOS 

183.9*" 

26.6 

26.8 

25.4 

25.2 

CANDELS EGS 

192.4'" 

26.6 

26.9 

25.5 

25.3 

Total 

896.3 






^ Photometry over a 31.8 and 34 arcmin^ area within the GS 
and GN fields is not available in the Bouwens et al. (2015) 
catalogs, due to these catalogs only including regions which 
have > 70% of the full depth available in the B 435 , VgOGi A 75 , 
zssOi ^I05i >1125, and bfieo bands. For these regions, we make 
use of the Skelton et al. (2014) catalogs to search for 2 > 7 
galaxies. 

^ The Bouwens et al. (2015) catalogs only cover the —450 
arcmin^ region from the CANDELS-UDS, COSMOS, and EGS 
fields where deep ACS and WFC3/IR data are available from 
CANDELS (75% of the area). In searching the CANDELS 
UDS, COSMOS, and EGS fields for 2 — 8 candidates, we make 
use of the Bouwens et al. (2015) catalogs where available and 
the Skelton et al. (2014) catalogs otherwise. 



photometric redshift 

Fig. 1.— Spitzer/IRAG [3.6] —[4.5] color vs. photometric redshift 
plot for young (—5 Myr) stellar population with very strong neb¬ 
ular emission lines (EW/fQ, = 1500A) and a flat continuum. Also 
ass umed are fixed flux ratios between em ission lines from Table 
1 of lAnders fc Fritze-v. AlvenslebenI II2003I ~) for 0.2 Zq metallicity, 
whilst assuming case B recombination for the Ha/H/3 flux ratio. 
The [3.6] — [4.5] color of galaxies is expected to become quite red 
at 2 > 7 due to the impact of the [OIII] line on the 4.5;rm band 
and no comparably bright nebular line in the 3.6/rm. 

the z > 7 LF. Finally, §6 includes a summary of our pa¬ 
per and a prospective. Throughout this paper, we refer 
to the HST F606W, F814W, F105W, F125W, F140W, 
and F160W bands as Feoe, .^ 814 , H 105 , J 125 , JHuo and 
H 160 , respectively, for simplicity. We also assume Hq = 
70 km/s/Mpc, = 0.3, and Da = O-'i"- A H magnitudes 
are in the AH system (|Oke fc Gunnlll983l) . 
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2. OBSERVATIONAL DATA SETS, PHOTOMETRY AND 
Z--5-8 SAMPLE 

2 . 1 . HST + ground-based data set and photometry 

The sample of z ^ 8 galaxies we identify in this paper 
is based on HST + ground-based observati ons that were 
acquired over 5 CANDELS a nd ERS fields (|Grogin et al.l 
120111 : HCoekemoer et al.l[2011l : Windhorst et al. 2011). 

The near-IR HST observations over the CANDELS 
fields range in depth from ^4 orbits over the ~130 
arcmin^ CANDELS DEEP components in COODS- 
North (CN) and COODS-South (GS) to 1 orbit depth 
over the ~550 arcmin^ CANDELS WIDE component in 
the GN, GS, UDS, COSMOS, and EGS fields. Over the 
GN and GS fields, the near-IR imaging observations are 
available in the Yios, J 125 , and Hiqq bands, while in the 
UDS, COSMOS, and EGS fields, the near-IR observa¬ 
tions are available in the J 125 and Hiqq bands. 

These fields also feature observations at optical wave¬ 
lengths with the HST ACS camera in the R 435 , 
beoe, * 775 ) -^814 and zgso bands for CANDELS-GN-I-GS 
(with 3-10-1- orbits per band), as well as Vgoe and 
7814 observations (^2-orbit depth) for the CANDELS- 
UDS-LCOSMOS-LEGS fields. 

In addition to the HST observations, these fields 
also have very deep ground-based observations from 
CFHT, Subaru Suprime-Cam, VLT HAWK-I, and 
VISTA/VIRCAM over the latter fields. Optical data are 
available in CANDELS-COSMOS field in the u, g, r, i, 
y and z bands as part of the CEHT legacy survey, and 
also in the B, g, V, r, i and z bands from Suby u ob¬ 
servations over the same field (|Capak et al.l 1200711 . The 
CANDELS-EGS field is observed in the same bands as 
the COSMOS field, as part of the CFHT legacy survey, 
whilst the CANDELS-UDS field is observed by Subaru 
as part of the Subary XMM-Newton Deep Field (SXDF) 
program (Furusawa et al. 2008). For extended sources, 
these optical observations reach similar or greater depths 
to the available HST data over these fields (i.e., 26 mag 
to 28 mag at 5cr in 1.2"-diameter apertures: see Bouwens 
et al. 2015) and allow us to exclude any potential lower 
redshift contaminants from our samples. 

Importantly, our ground-based observations also in¬ 
clude moderately-deep (~26 mag at 5cr [1.2"-diameter 
apertures]) Y-band observations which we use to con¬ 
strain the nature of our selected z > 7 can¬ 

didates (where HST observations are unavailable). 
These observations are available over the CANDELS- 
UDS-I-COSMOS fi elds through HAWK -I and VISTA as 
part of the HUGS (iForhana et al.l[20lll and UltraVISTA 
(|McCracken et al.ll2012f) programs, respectively. A more 
detailed description of the observations we utilize in con¬ 
structing our source catalogs, as well as o ur procedure for 
constr ucting these catalogs is provided in Bouwens_^t_al. 
(1201^ (see Table 1, Figure 2, and §3 from iBouwens et 

[ml . 

HST photometry was performed running t he Source 
Extractor software ([Bertin fc Arnouti 1199611 in dual¬ 
image mode, taking the detection images to be the 
square root of the image (iSzalav et al.lll999ll and PSF- 
matching the observations to the Higo-band PSF. The 
colors and total magnitudes were measured with Kron- 
like (1980) apertures and Kron factors of 1.6 and 2.5 
respectively. 


I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I 

[ 3 . 6 ]-[ 4 . 5 ] > 0,5 
selection 



6.0 6.5 7.0 7.5 8.0 8.5 9.0 9.5 

photometric redshift 

Fig. 2.— The photometric redshift distribution of the 2 = 7- 
9 [3.6] — [4.5] > 0.5 IRAC-selected control sample (including 15 
sources) we used to validate our HST-hlRAC selection technique 
(§3.2). The control sample of IRAC-red optical dropouts was 
identified exclusively from those fields with deep T-band data 
(i.e., CANDELS GOODS-North, GOODS-South, UDS, COSMOS, 
ERS) and consider sources that are fainter than we focus on here 
for our primary selection (Table [^. Redshift estimates were made, 
based on their observed HST - 1 - ground-based photometry. No 
consideration of the Spitzer/IRAC fluxes is made in deriving the 
photometric redshift presented here (to ensure that the two red¬ 
shift measures are entirely independent). These results strongly 
suggest that one can use the Spitzer/IRAC [3.6] —[4.5] color to re¬ 
liably distinguish 2 > 7 galaxies from z < 7 galaxies (especially 
in the present case where one makes exclusive use of those sources 
with relatively extreme [3.6] — [4.5] colors). 

Photometry on sources in the ground-based data is per¬ 
formed after the contamination from foreground sources 
is removed, using an automated clea ning procedure 
(jLabbe et al.ll201(M iLabbe et al.ll2010bh . The positions 
and two-dimensional spatial profiles of the foreground 
sources are assumed to match that seen in the high- 
spatial resolution HST images, after PSF-matching to 
the ground-based observations. The total flux in each 
source is then varied to obtain a good match to the light 
in the ground-based images. Light from the foreground 
sources is subsequently subtracted from the images, be¬ 
fore doing photometry on the sources of interest. Flux 
measurements for individual sources are then performed 
in 1.2"-diameter circular apertures due to the objects 
being inherently unresolved in the ground-based obser¬ 
vations. These flux measurements are then corrected to 
total, based on the model flux profiles computed for in¬ 
dividual sources based on the observed PSFs. The pro¬ 
cedure we employ here to derive fluxes is very similar to 
that employed in Skelton et al. (2014: see also Galametz 
et al. 2013 and Guo et al. 2013 who have also adopted a 
similar procedure for their ground-based photometry). 

2.2. Spitzer/IRAC Data Set and Photometry 

The detailed information we have on z ^ 6-9 galaxy 
candidates over the GANDELS fields from HST is nicely 
complemented in the mid-IR by the Spitz er Extended 
Deep Survey (SEDS, PI: Fazio) program (jAshbv et al.l 
l2013f) . which ranges in depth from 12 hours to >100 
hours per pointing, though 12 hours is the typical ex- 
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posure time. The SEDS program provides us with flux 
information at 3.6/rm and 4.5/rm, which can be useful for 
probing 2 ~ 6-9 galaxies in the rest-frame optical, quan¬ 
tifying the flux in various nebular emission lines, and 
estimating the redshift. 

Over the GOODS-North and GOODS-South fields, we 
make use of Spitzer/IRAC reductions, which include es¬ 
sentially all the Spitzer/IRAC observations obtained to 
the present (Labbe et al. 2015: but see also Ashby et 
al. 2015), with 50-200 hours of observations per pixel in 
both bands (and typically ^100 hours). 

Our procedure for performing photometry on the 
IRAC data is essentially identical to that used on the 
ground-based observations, except that we utilize 2"- 
diameter circular apertures for measuring fluxes. These 
fluxes are then corrected to total based on the model 
profile of the individual sources -I- the PSF. Depending 
on the size of the source, these corrections range from 
^2.2x to 2.4x. 

The median 5cr depths of these Spitzer/IRAC obser¬ 
vations for a ^26-mag source is 25.5 mag in the 3.6^m 
band and 25.3 mag in the 4.5/rm band. 


erwise, we will rely on the Skelton et al. (2014) catalogs 
and photometry. 

We then apply color criterion to identify a base sample 
of Lyman-break galaxies at z ^ 6.3-9.0. In particular, 
over the CANDELS-UDS, COSMOS, and ECS fields, we 
use a 

(.^814 — Jl 2 b > 2.2) A (J 125 ~ Hiqq < 0.5)A 
(.^814 — Jl2b > 2(Ji25 — Hiqo) + 2.2) 

criterion. Over the CANDELS-GN and GS fields, we 
require that sources satisfy one of the two color criteria 
defined by Eq. [2] or Eq. [H 

(^850 ~ ^105 > 0-7) A (Ji25 — Hiqo < 0.45)A 
(^850 ~ ^105 > 0.8(Ji25 — Hiqo) + 0.7)A (2) 

((^814 — Jl25 > 1-0) V {SN{Isii) < 1.5)) 

(^105 ~ Ji25 > 0.45) A (Ji25 — Hiqo < 0.5) 

A (^105 ~ Ji25 > 0.75( Ji25 ~ Hiqo) + 0.525) 


3. SAMPLE SELECTION 


3.1. [3.6]-[4-5] IRAC color vs. redshift and HST 
detections 


Many recent st u dies (e.g., iSchaerer fc de Barrod 


Stark et ^ 


Shim.et_^n 200; jSmjt_et_aij_ |2014 ;_ ^bbe et al.l 


2009; 


2013; 


2013t Ide Barros et al.l 1201411 have presented 


convincing evidence to support the presence of strong 
nebular line contamination in photometric filters, partic¬ 
ularly for the Spitzer/IRAC [3.6] and [4.5] bands. The 
observed [3.6] —[4.5] IRAC color of galaxy candidates ap¬ 
pear to be strongly impacted by the presence of these 
lines at different redshifts, in particular those of Ha and 
[OIII]. Figure [U provides an illustration of the expected 
dependence of the Spitzer/IRAC [3.6] — [4.5] color, as a 
function of redshift, assuming an [OIII]-fH// EW (rest- 


frame) of ^2250A, which is at the hi gh end of what has 
been estimated for galaxies at z 7 (jhabbe et al.ll201,^ 
iSmit et ^I20l4 iSmit et ^120151) . 

The significant change in the [3.6] — [4.5] color of galax¬ 
ies from z ^ 6-7 to z > 7 suggests this might be a promis¬ 
ing way of segregating sources by redshift and in particu¬ 
lar to identify galaxies at z > 7. Such information would 
be especially useful for search fields like CANDELS ECS, 
which lack deep observations in E-band at ^l.l/rm to es¬ 
timate the redshifts directly from the position of the Ly¬ 
man break. Smit et al. (2015) have shown that selecting 
sources with blue [3.6] — [4.5] colors can effectively single- 
out sources at z 6.6-6.9 over all CANDELS fields, even 
in the absence of E-band coverage. 

Here we attempt to exploit this strong dependence of 
the [3.6] — [4.5] color on redshift to identify some of the 
brightest z > 7 galaxies over the CANDELS fields. In 
performing this selection, we start with the source cata¬ 
logs derived by Bouwens et al. (2015) and Skelton et al. 
(2014) over a ^900 arcmin^ region from the five CAN- 
DELS fields. In general, we will rely on the source cat¬ 
alogs from Bouwens et al. (2015) where they exist (cov¬ 
ering a 750 arcmin^ area or ^83% of CANDELS) 0 0th- 


® Bouwens et al. (2015) only considered those regions in CAN- 


These color criteria are essentially identical to those 
from Bouwens et al. (2015) but allow for J 125 — Hiqo 
colors as red as 0.5 mag to match up with the color cri¬ 
teria of Oesch et al. (2014) and Bouwens et al. (2015) in 
searching for z > 8.5 galaxies (i.e., J 125 —Llieo > 0.5). In 
so doing, it was our goal to maximize the completeness 
of our selection for bright z = 7-9 galaxies within the 
CANDELS program[3 

These color criteria are motivated in Figure 3 of 
Bouwens et al. (2015) and result in a very similar redshift 
segregation as one achieves using photometric redshifts. 

We require that sources have [3.6] — [4.5] colors redder 
than 0.5 mag (see Figure [T|). This color criterion was 
chosen (1) so as to require slightly redder colors than 
was the average color measured by Labbe et al. (2013) 
for their faint z ~ 8 sample from the HUDF (i.e., ~0.4) 
and (2) such that sources would not easily satisfy the 
criterion simply due to noise (requiring >2a deviations 
for the typical source). To be certain that the IRAC 
colors we measured are robust, we exclude any sources 
where the subtracted flux from neighboring sources ex¬ 
ceeds 65% of the original flux in a 2"-diameter aperture 
(before subtraction). 

To ensure our selection is free of z < 7 galaxies, we re¬ 
quired that sources show no statistically significant flux 
at optical wavelengths. Sources that show at least a 
1.5(7 detection in terms of the inverse-variance-weighted 
mean Vqqq and /814 flux with HST were excluded. In 
addition, we also excluded sources detected at > 2.5(7 
in the deep optical imaging observations available over 
each field from the ground. We adopted a slightly less 
stringent threshold for detections in the ground-based 
observations, due to the impact of neighboring sources 
on the overall noise properties. 

Finally, we consider potential contamination by low- 
mass stars, particularly later T and E dwarfs (T4 and 

DELS where deep optical and near-IR observations are available 
from the CANDELS observations. 

We remark that any contaminants in a particularly bright 
selection would be generally easy to identify given the depth of the 
HST, Spitzer, and supporting ground-based observations. 
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later), where the [3.6] — [4.5] color can become quite a 
bit redder than 0.5 mag (Kirkpatrick et al. 2011; Wilkins 
et al. 2014). To exclude such sources from our samples, 
both the spatial information we had on each source from 
the SExtractor stellarity parameter and total SED in¬ 
formation were considered. Sources with measured stel- 
larities >0.9 were identified as probable stars (where 0 
and 1 corresponds to extended and point-like sources, 
respectively), as were sources with measured stellarity 
parameters >0.5 if the flux information we had avail¬ 
able for sources was significantly better fit (Ay^ > 2 ) 
with a low-mass stellar model from the SpecX prism 
library (Burgasser et al. 2004) than with the best-fit 
galaxy SED model, as derived by the Easy and Ac- 
curate Zphot from Yale (EAZY; iBrammer et al.l [2008fl 
software. Our SED fits with EAZY considered both 
the standard SED templates from EAZY and SED tem¬ 
plates fro m the Galaxy Evolu tionary Synthesis Models 
(GALEV; iKotulla et all 1200911. Nebular emission lines, 
as described by lAnders fc Fritze-v. AlvenslebenI (1200311 . 
were added to the GALEV SED template models assum¬ 
ing a 0 . 2^0 metallicity. 

No sources were removed from our selection as proba¬ 
ble low-mass stars. The procedure we use here to exclude 
low-mas s stars from our select ion is identical to that uti¬ 
lized bv iBouwens et al.l (1201511 . 

3.2. Validation of Selection Technique 

Before applying the selection criteria from §3.1 to the 
^900-arcmin^ GANDELS -I- ERS search fields, it is use¬ 
ful to first test these criteria on those data sets which 
feature deep x and Y-band observations. The availabil¬ 
ity of observations at these wavelengths, together with 
observations at both redder and bluer wavelengths with 
HST, allows for very accurate estimates of the redshifts 
for individual sources. There are five data sets that pos¬ 
sess these observations: (1) GANDELS GOODS-S, (2) 
GANDELS GOODS-N, (3) ERS, (4) GANDELS UDS, 
and (5) GANDELS COSMOS field. The first three fea¬ 
ture these observations with HST and the latter two us¬ 
ing ground-based telescopes. 

We apply selection criteria from the previous section 
to a iLi 6 o-band limiting magnitude of 26.7 mag for the 
first three fields and 26.5 mag for the latter two. Our 
decision to use these depths is partially guided by the 
sensivity of the Spitzer/IRAC data over these fields. 

Applying the selection criteria from the previous 
section to the GANDELS GN-I-GS and ERS fields 
{Hi&o,ab < 26.7), we find 7 sources that satisfy our se¬ 
lection criteria. For each of these sources, we estimate 
photometric redshifts with EAZY. In fitting to the ob¬ 
served photometry, we used the same standard EAZY 
SED templates as we described in the previous section. 

We also applied the above selection criteria to the 
CANDELS-UDS and CANDELS-COSMOS fields, where 
it is also possible to estimate photometric redshifts, mak¬ 
ing use of the available HST observations and ground- 
based optical and near-IR Y and K band observations. 
Eight sources satisfy these criteria. 

All 15 of the sources selected using the criteria from 
the previous section are presented in Figure [5] and fall 
between z = 7.0 and z = 8.3, which is the expected 
range if a high-EW [OHI]-|-H/3 line is responsible for red 
[3.6] — [4.5] colors in these galaxies. This suggests that 


the criteria we propose in the previous section can be ef¬ 
fective in identifying a fraction of z > 7 galaxies that are 
present in fields with deep HST-(-Spitzer observations. 
The individual coordinates, colors, and estimated red¬ 
shifts for individual sources from this validation sample 
can be found in Table [5| located in Appendix B. 

In recommending the use of the IRAG photometry to 
subdivide z 6-9 samples by redshift, we should empha¬ 
size that the most robust results will be obtained making 
use of only those sources with the smallest confusion cor¬ 
rections. While we took care in the selection of both our 
primary sample (and the sample we used to validate the 
technique) to avoid such souces, such sources were not 
excluded in making Figure 1 of Smit et al. (2015: re¬ 
sulting in a few z > 7 sources with anomalously blue 
Spitzer/IRAG colors). Despite this issue with Smit et 
al. (2015) Figure 1, we emphasize that this is neverthe¬ 
less not a major concern for sources in their z = 6 . 6 - 6 .9 
sample. Only 2 of the 15 sources in the latter sample 
were subject to a ^3x correction for flux from neighbor¬ 
ing sources and those 2 sources (GSD-2504846559 and 
EGS-1350184593) are flagged as less reliable. 

3.3. Search Results for Bright Hiqq ab < 25.5 Galaxies 

Here we focus on the identihcation of only the brightest 
Hi60,ab < 25.5 z > 7 galaxies using our Spitzer/IRAG 
color criteria. This is to keep the current selection small 
and to focus on sources whose surface density was par¬ 
ticularly poorly defined from previous work. Prior to 
this work, the only study which identified such bright 
z ^ 8 sources was Bouwens et al. (2015). Focusing on 
the brightest sources is also valuable, since it allows us 
to obtain very precise constraints on SED shapes and its 
Spitzer/IRAG colors of the sources, as well as providing 
opportunities for follow-up spectroscopy (see §4. 2). 

Applying the selection criteria described in H3.ll on 
the GANDELS-GS, CANDELS-GN, GANDELS-UDS, 
CANDELS-GOSMOS and GANDELS-EGS fields, we 
identify a total of 4 especially bright (Hieo.AS < 25.5) 
candidate z >7 galaxies. 

Our 4 candidate z >7 galaxies are presented in Table 
[2] and in Figure[3l We see from Figure[3|that each candi¬ 
date is clearly visible in the HST Hieo and J 125 filters, as 
well as the IRAG 3.6 ^m and 4.5 /rm bands. Of course, 
no signihcant detection is evident in the HST Feoe and 
1814 bands for these sources. This would suggest that 
these sources show a break in their spectrum somewhere 
between 0.9/rm and 1.2^m and therefore have redshifts 
between z ^ 6 and z ~ 8.5. [We discuss the impact of 
information from Y -band observations available over 3 of 
the 4 candidates in §4.1.] 

3 of these 4 bright sources are found in the GANDELS- 
EGS held. Sources from this held were not included in 
our earlier attempt to validate the present selection tech¬ 
nique (§3.2), so only one of these new sources is in com¬ 
mon with the 15 sources just discussed. 

To derive constraints on the redshift of each bright 
source, we again made use of EAZY. The photome¬ 
try provided to EAZY included huxes from HST hlters, 
IRAG 3.6 /rm and 4.5 ^m hlters and ground based tele¬ 
scopes. Using EAZY allows us to generate a best-ht SED 
of each galaxy candidate as well as its redshift likelihood 
distribution (P(z)), which we present in Figured with 
the observed galaxy hux points overplotted. From the 
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TABLE 2 

A COMPLETE LIST OF THE RESULTING Z>7 SOURCES IDENTIFIED AFTER APPLYING OUR SELECTION CRITERIA. 


ID 

R.A. 

Dec 


[3.6]-[4.5] 

^phot^ 

Dos — 

References* 

COSY-0237620370 

EGS-zs 8 -l 

EGS-ZS8-2 

EGSY-2008532660 

10:00:23.76 

14:20:34.89 

14:20:12.09 

14:20:08.50 

02:20:37.00 

53:00:15.35 

53:00:26.97 

52:53:26.60 

25.06±0.06 

25.03±0.05 

25.12±0.05 

25.26±0.09 

1.03±0.15 

0.53±0.09 

0.96±0.17 

0.76±0.14 

7 -| 4 _|_ 0.12 

7.92±3:| 

-0.13±0.66 

1.00±0.60 

0.66±0.37 

[1],[2],[3] 

[3], [4] 

[3] 


* References: [1] Tilvi et al. 2013, [2] Bowler et al. 2014, [3] Bouwens et al. 2015, [4] Oesch et al. 2015b 
^ The apparent magnitude of each source in the /fieo band. 

The photometric redshift estimated by EAZY, including flux measurements in the Y band. The uncertainties 
quoted here correspond to Icr. 

The Y — J color for each source. The COSMOS candidate uses ground based data whilst the EGS candidates use 
Yios and J 125 filters (where available). 



Fig. 3.— HST/ACS V606^814, HST/WFC3 Yios Ji257fi60, and Spitzer/IRAC 3.6/rm+4.5/rm postage stamp images (4" x4”) of the 3 z >7 
candidates identified over the 5 CANDELS fields. On the Spitzer/IRAC images, flux from neighbouring sources has been removed. Y-band 
observations at 1.05;rm are also available for COSY-0237620370 from ground-based programs (ZFOURGE [Tilvi et al. 2013], UltraVISTA 
[Bowler et al. 2014]). 


SED plots, we observe a near-flat rest-frame optical con¬ 
tinuum, as well as emission lines dominating at the loca¬ 
tion of high flux points, highlighting the contribution of 
strong nebular emission lines to the instrument filters. 

One of our z > 7 candidates, i.e., EGS-zs8-2, is suf¬ 
ficiently compact, as can be seen from Eigure |31 that 
we considered the possibility that it may correspond to 
a star. To test this possibility, we compare its SED to 
all the stellar SEDs in the SpecX prism library and find 
the best-fitting stellar SED. The goodness-of-fit for 
the stellar SED is an order of magnitude greater than 
the galaxy SED. In addition, the SExtractor stellarity 
we measure for EGS-zs8-2 in the J 125 and i/ieo bands is 
0.60 and 0.33 (where 0 and 1 correspond to an extended 
and point source, respectively), which significantly favors 
EGS-zs8-2 corresponding an extended source. Bouwens 


et al. (2015) ran an extensive number of end-to-end simu¬ 
lations to test the possibility that point-like sources could 
scatter to such low measured stellarities. Stellarities of 
~0.60 are only found for Hiqq ab 25-mag point-like 
sources in <2% of the simulations that Bouwens et al. 
(2015) run. Therefore, both because of the spatial and 
spectral information, we can be confident that the EGS- 
zs8-2 candidate is a z > 7 galaxy and not a low-mass star. 
However, as we show in §4.2, perhaps the most convinc¬ 
ing piece of evidence for this source corresponding to a 
z > 7 galaxy is our discovery of a plausible 4.7a Lya line 
in the spectrum of this source at 1.031/Lm (Eigure [5]). 

A second candidate from our selection, COSY- 
0237620370, is also very compact and could poten¬ 
tially also correspond to a low-mass star. However, 
like EGS-zs8-2, the photometry of the source is bet- 
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A (jLim) z 


, SED models (blue line) to the observed HST + Spitzer/IRAC + ground-based photometry {red points and error 
r bright (-H^i 60 ,AS < 25.5) z > 7 galaxies selected using our IRAC-red selection criteria ([3.6] — [4.5] > 0.5). Also 
the redshift estimate for the best-fit model SED provided by EAZY. Right: Redshift likelihood distributions P{z) 
^ z > 7 galaxies, as derived by EAZY. The impact of the Spitzer/IRAC photometry on the redshift likelihood 
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Fig. 5. — Redshift likelihood P{z) constraints for the respec¬ 
tive galaxy candidates presented in Figure [4| shown considering 
the impact of various subsets of the photometry (i.e., HST only, 
HSTH-IRAC, HST-|-IRAC-|-ground-based observations (no Y-band 
and Y-band included). It is clear that the addition of both the 
Spitzer/IRAC observations and the Y-band observations result in 
a much tighter distribution and allows for a much more accurate 
estimation of the photometric redshift. 


ter fit with a galaxy SED than a stellar SED (with a 
X^istar) — (galaxy) = 17.2) and the source shows ev¬ 
idence for spatial extension, with a measured stellarity 
of 0.81 and 0.34 in the J 125 and Tfieo bands, respec¬ 
tively. Stellarities even as high as 0.81 are only recov¬ 
ered in ^5% of the end-to-end simulations Bouwens et al. 
(2015) run at iJieo.AB-band magnitudes of ^25.0. Ear¬ 
lier, Tilvi et al. (2013) obtained the same conclusion re¬ 
garding this source based on medium-band observations 
over this candidate from the ZFOURGE program, where 
consistent fluxes are found in the near-infrared medium 
bands strongly arguing against this source correspond¬ 
ing to a low-mass star. Bowler et al. (2014) also con¬ 
clude this source is extended and not a low-mass star, 
based on its spatial profile (see Figure 6 from Bowler et 
al. 2014) and based on its observed photometry where 
y^(star) — x^ (galaxy) = 13.0. 

Flux information from HST, Spitzer/IRAC, and 
ground-based observations all have value in constrain¬ 
ing the redshifts of the candidate z >7 galaxies we have 
identified in the present probe. While the HST flux in¬ 
formation we have available for all three candidates in 
the VeoQ IS 14 J 125 JHno H 160 bands only allow us to place 
them in the redshift interval z ^ 6 .5-9.0 (blue line in 
Figure 0 , we can obtain improved constraints on the 
redshifts of the candidates incorporating the flux infor¬ 
mation from Spitzer/IRAC and from deep ground-based 
observations. Each of these three candidates appears to 
have redshifts robustly between z ^ 7.0 and z ^ 8.6 (red 
and black lines in Figure0. In addition, as we discuss in 
§4.1 and show in Figure[Sl the availability of the Y /Yios- 
band observations allows us to significantly improve our 
redshift constraints on all three candidates. 

We used the Bouwens et al. (2015) catalogs to search 
83% of the total area of the CANDELS fields and the 
Skelton et al. (2014) photometric catalogs otherwise (in 
those regions over the WFC3/IR CANDELS fields which 
lack the deep HST/ACS data). As a check on the search 
results we obtained with the Bouwens et al. (2015) cata¬ 
logs, we applied the same selection criteria to the Skelton 
et al. (2014) catalogs. Encouragingly, we identified 75% 
of our sample, with only one candidate missing due to its 
having a [3.6] — [4.5] color of 0.47 mag. For all 4 candi¬ 
dates from our primary sample, we find that our derived 
[3.6] —[4.5] colors are almost identical to those quoted by 
Skelton et al. (2014), agreeing to <0.1 mag (and typically 
A[3.6]-[4.5] of 0.05 mag). 

We also identified 1 additional bright (iJieo.AS < 25.5) 
z > 7 candidate in the CANDELS-ECS field not iden¬ 
tified in our primary search (see Appendix A). It seems 
clear examining its photometry that this source is ex¬ 
tremely likely to be at z ~ 7-9 (and indeed it appears in 
the Bouwens et al. 2015 z ^ 8 sample). However, since 
its measured [3.6] — [4.5] color is 0.22±0.06 mag in our 
photometric catalog (0.3 mag bluer than in the Skelton 
et al. 2014 catalog), we did not include it in our primary 
sample. We remark that photometry for this source was 
more challenging due to its being located close to a bright 
neighbor and its being a two-component source. 


3.4. Possible Evidence for Lensing Amplification of 
Selected z > 7 Sources 
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Fig. 6 .— Image of the area (9.6” X 7.2”) surrounding our new 
z ~ 8.6 candidate EGSY-2008532660. EGSY-2008532660 lies very 
close (< 3”) to two bright, apparently massive foreground galax¬ 
ies (§3.4). Based on the position of the foreground sources and 
their inferred masses, we estimate that EGSY-2008532660 is likely 
magnified by a modest factor, i.e., ~ 1.8x. 

For very high redshift sources {z » 6), it is ex¬ 
pected that the sources with the brightest apparent mag¬ 
nitudes will benefit from gravitational lensing (Wyithe et 
al. 2011; Barone-Nugent et al. 2015; Mason et al. 2015; 
Fialkov & Loeb 2015), and indeed it is found that a 
small fraction of the brightest galaxies identified over the 
CANDELS program are consistent with being boosted by 
gravitational lensing (Barone-Nugent et al. 2015). 

To investigate whether any of the bright z >7 galaxies 
identified in our search might be gravitationally lensed, 
we considered all sources within 5" of our candidates in 
the Skelton et al. (2014) catalogs and used the estimated 
redshifts, stellar masses, and sizes from these catalogs to 
derive Einstein radii for the foreground sources assuming 
a single isothermal sphere model. We then calculated the 
degree to which our bright z >7 galaxy candidates might 
be magnified by the foreground sources. 

In only one case was the expected magnification level 
>10% and this was for our z ~ 8.6 candidate EGSY- 
2008532660. In this case, we identified two foreground 
galaxies which could significantly magnify this candidate 
(Figure HI). The first was a Mq mass, z ^ 1.4 

galaxy (14:20:08.81, 52:53:27.2) with a separation of 2.8" 
from our z ~ 8.6 candidate. The second was a 10^^'^ 
Mq mass, z ~ 3.1 galaxy (14:20:08.37, 52:53:29.1) with 
a separation of 2.7" from our z ~ 8.6 candidate. Using 
the measured size of the two sources, we derive u ~ 170 
km/s and a ~ 370 km/s for the velocity dispersion. We 
checked and these velocity dispersions are fairly similar 
to what fitting formula in Mason et al. (2015) yield (i.e., 
using the relation in their Table 1 and applying iJieo- 
band or IRAC 3.6/im apparent magnitudes depending 
on whether we are considering the z~1.4orz^3.1 
source). 

Based on the observed separation of this source from 
our z ^ 8.6 galaxy, we estimate a lensing magnification 
of 20% and a factor of 1.8 from the former and latter 
foreground sources. In computing these magnification 
factors, we assume that the mass profile of galaxies is an 
isothermal sphere and taking the magnification factor to 



z 


Fig. 7.— Redshift likelihood distributions P(z) for the 4 
t7i60 AB < 25.5 sources in our 2 > 6, IRAG ultra-red selection. 
These likelihood distributions include flux constraints from the 
ground-based and HST Y-band observations. Our three highest 
redshift sources have = 7.6 ± 0.3, = 7.9 ± 0.4, and 

^phot — 8.6itQ 4. 

be 1/(1 — OeIO) where 9 is separation from the neighbor¬ 
ing sources and 9e is the Einstein radius. Looking at the 
morphology of EGSY-2008532660, we see no clear evi¬ 
dence to suggest that the galaxy is highly magnified and 
there is no obvious counterimage. However, we clearly 
cannot rule out smaller lensing amplification factors, par¬ 
ticularly if the intrinsic size of the source is small. As the 
inferred stellar or halo masses for the neighboring galax¬ 
ies is not precisely known, this translates into a modest 
uncertainty into the actual luminosity of this source (as 
much as 0.3 dex). Given this fact, we consider it safest 
for us to exclude it from analyses of the UV LF. 

4. VALIDATION OF OUR Z ~ 8 SELEGTION 

Here we attempt to determine the nature of 
the z > 7 candidates we selected using the 
HST+Spitzer/IRAG+ground-based observations using 
some Y-band observations that became available over 
a few of our candidates and using the results of some 
follow-up spectroscopy that we performed (first reported 
in Oesch et al. 2015b). 

4.1. Y-band Photometric Observations 

Deep observations at 1.05/im are particularly useful in 
ascertaining the nature of these candidates and also their 
redshift, due to the Y-band photometry providing con¬ 
straints on the position of the Lyman break as it redshifts 
from 1.2/rm to 0.9/rm. 

Deep observations at 1.05/im are available for 3 of the 
4 z > 7 candidates that we selected as part of our 
Hi 60 ,ab < 25.5 sample. Y-band observations of the 
GOSY-0237620370 candidate are available from the 3- 
year UltraVlSTA observations (McCracken et al. 2012), 
while HST Yios-band observations are available over 2 
other candidates in our selection as a result of some re¬ 
cent observations from the z9-CANDELS follow-up pro- 
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gram (Bouwens 2014; Bouwens et al. 2016)E3 

We make these estimates in an identical way to what 
we did previously. Our redshift constraints, including 
the Wband, are presented in Table [2] Furthermore, we 
present the HST Y 105 filter images in Figure [31 where 
we observe a clear detection in the Y 105 filters for EGS- 
zs 8 -l and EGS-zs8-2 and no detection in the Feoe or 
/814 filters, indicating a 2~7 Lyman dropout. For EGS- 
zs 8 -l, however, we observe little to no detection in the 
Y 105 filter but a clear detection in the J 125 filter which 
indicates this galaxy is observed at 2 ~ 8 . 

Figure [Hand [7] presents the redshift likelihood distribu¬ 
tions on our z >1 candidates, incorporating the E-band 
observations from UltraVlSTA and HST. It is evident 
from Figure [5] that the E-band data greatly improves 
our constraints on the redshift of the individual candi¬ 
dates in our selection. Together with the results in §3.2 
and Figure 121 these results largely validate our selection 
technique. 

4.2. Keck/MOSFIRE Spectroscopic Follow-up 
4.2.1. Observations and Reduction 

In addition to using photometric data in the Y-band 
to validate our method, we also tested this method by 
obtaining deep near-IR spectroscopy on 2 sources from 
the current selection. Oesch et al. (2015b) already pro¬ 
vided a first description of the observational set-up we 
utilized for half of our targets, so we keep the current 
discussion short. A total of 4 hours of good Y-band 
spectroscopy were obtained in the CANDELS-EGS field 
with the Multi-Object Spectrometer for Infra-Red Ex¬ 
ploration (MOSFIRE: McLean et al. 2012) instrument 
on the Keck I te lescope. Two masks (see Fig. 2 of 
lOesch et al.l 1201511 were utilized and our spectra were 
taken with 180 s exposures at a spectral resolution of 
R=3500 and R=2850 (for a 0.7” and 0.9” slit respec¬ 
tively) over 3 nights (April 18, April 23, April 25, 2014 
- although due to poor weather conditions, April 18 was 
effectively lost), with the aim of searching for Lya emis¬ 
sion in EGS-zs 8 -l and EGS-zs8-2. Each mask contains a 
slitlet placed on a star, which we use for monitoring the 
sky transparency and observing conditions of each expo¬ 
sure. These observations were reduced using a modified 
version of the DRP MOSFIRE reduction code pipeline 
(for details see Oesch et al. 2015b). The spectra com¬ 
plement the photometric data sets for these two galaxies 
and allow us to confirm their redshifts. 

4.2.2. Lya Emission Lines 

The observations carried out with Keck/MOSFIRE re¬ 
vealed candidate Lya emission lines in the spectra of 
both EGS-zs 8 -l and EGS-zs8-2. The detection of a Lya 
line for EGS-zs 8 -l appears to be robust (a 6.1cr detec¬ 
tion with a line flux of fhyct = 1-7 ± 0.3 x 10“^^ erg s“^) 
and places that source at ZLya=7.7302±0.0006, as first 
reported by Oesch et al. (2015b)13 

The purpose of the z9-CANDELS program was to determine 
the nature of high-probability but uncertain candidate 2 9-10 

galaxies over the CANDELS-UDS, COSMOS, and EOS fields. In 
some cases, bright candidate 2 ~ 8 galaxies were located nearby 
bright 2 ~ 9-10 candidates and could be readily observed in the 
same pointings. 

The flux uncertainties that we derive for this candidate and 
EGS-zs8-2 is almost an order-of-magnitude larger than found in ob- 


The ID and 2D spectra for our other targeted z > 7 
candidate EGS-zs8-2 is presented in Figure [ 8 ] (see also 
Figure 2 in Oesch et al. 2015b for spectra of the confirmed 
z = 7.7302 ± 0.0006 candidate). Using a simple gaussian 
to determine the central wavelength of the observed line 
at 1.031/im (and ignoring asymmetry and other effects 
due to skylines surrounding this candidate Lya line), we 
determine the spectroscopic redshift for the source to be 
ZLya = 7.4770 ± 0.0008, with a detection significance of 
4.7 ct for the line and a line flux of fLya = 1-6 ± 0.3 x 
10“^^ erg s“^ cm“^. While this line is only detected at 
4.7cr significance, its reality appears to be supported by 
subsequent near-infrared spectroscopy obtained on this 
source from independent observing efforts (D. Stark et 
al. 2016, in prep). 

In addition to the Lya-emission lines reported by 
Oesch et al. (2015b) and this work, Zitrin et al. (2015) 
report the detection of a 7.5cr Lya line for our EGSY- 
2008532660 candidate in new Keck/MOSFIRE obser¬ 
vations (June 10-11, 2015). This redshift measure¬ 
ment sets a new high-redshift distance record for galax¬ 
ies with spectroscopic confirmation. Our photometric 
selection therefore contains 3 of the 4 most distant, 
spectroscopically-confirmed galaxies to date. 

The ZLya = 7.730, ZLya = 7.477 and ZLya = 
8.683 redshifts for EGS-zs 8 -l, EGS-zs8-2 and EGSY- 
2008532660, respectively, are in excellent agreement with 
the photometric redshifts derived for these galaxies using 
HST-|-IRAC-|-Ground based observations and our color 
criteria. The absolute magnitude and redshifts of EGS- 
zs 8 -l, EGS-zs8-2, and EGSY-2008532660 are presented 
in the top panel of Figure [9] in relation to other z > 6.5 
galaxies with clear redshift determinations from Lya. 

The current specroscopy provides us with considerable 
reassurance that our proposed color technique is an ef¬ 
fective method to search for bright, z > 7 galaxies 

5. COMPARISON WITH PREVIOUS WORK 

Three of our four candidates were already identified 
as part of previous work. Tilvi et al. (2013) identified 
GOSY-0237620370 as a z ^ 7 galaxy by applying Lyman- 
break-like criteria to the deep medium-band ZFOURGE 
data and estimate a redshift of 7.16lg^g. This source was 
also identified by Bowler et al. (2014) as a z ^ 7 galaxy 
(211127 in the Bowler et al. 2014 catalog) using the deep 
near-IR observations from the Ultra-VISTA program and 
derived a photometric redshift of 7.03lo;}j for the source 
(or 7.20 if the source exhibits prominent Lya emission), 
similar to what we find here. Tilvi et al. (2013) derive 

servations of similar 2 > 7 (e.g., Finkelstein et al. 2013). This is in 
part due to the significantly poorer seeing conditions we were sub¬ 
ject for the observations (1.00^^ FWHM instead 0.65^^ for Finkel¬ 
stein et al. 2013). Another potentially significant contributing fac¬ 
tor is our relatively conservative account of the uncertainties in 
the line flux measurements, including uncertainties that arise from 
the sky subtraction. The uncertainties we derive are consistent 
with typical values reported by the MOSDEF program (Kriek et 
al. 2015). 

Interestingly enough, D. Stark et al. (2016, in prep) have 
also spectroscopically confirmed that the fourth source (COSY- 
0237620370) from our sample lies at z=7.15. As such, Lyo emis¬ 
sion has been found in all 4 galaxies that make up our selection. 
Our entire sample has therefore been spectroscopically confirmed 
to lie in the redshift range 2 = 7.1-9.1, with the spectroscopic red¬ 
shifts being in excellent agreement with our derived photometric 
redshifts. 
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Fig. 8 .— Keck/MOSFIRE spectra of EGS-zs8-2. The 2D spectrum after a 2x2 binning is presented in the upper panel, while the extracted 
ID spectrum is shown in the lower panel. The gray shaded area represents the Icr flux uncertainty, the red line shows the best-fit Gaussian. 
A candidate Lya line (detected at 4.7cr significance) is apparent at 1.031/im between two sky lines. Using a simple gaussian to model the 
shape and position of this line suggests a redshift of z = 7.4770 zb 0.0008 for this source (see also D. Stark et al. 2016, in prep). The other 
z > 7 candidate here targeted with spectroscopy also shows a prominent Lya line, with a measured redshift of 7.7302 ± 0.0006 (Oesch et 
al. 2015b). 


a [3.6] — [4.5] color of 1.96 ± 0.54 mag, while Bowler et 
al. (2014) find 0.7±0.3 mag, both of which are broadly 
consistent with what we find here. 

Bouwens et al. (2015) identified 3 of the 4 sources as 
part of their search for z ^ 7-8 galaxies over the five 
CANDELS fields and segregated the sources into dif¬ 
ferent redshift bins using the photometric redshift esti¬ 
mates. The full HST -|- Subaru Suprime-Cam BgVriz 
+ CFHT Megacam ugriyz + UltraVISTA YJHKs 
photometry was used to estimate these redshifts for 
the candidate in the COSMOS field. Meanwhile, the 
HST -I- CFHT Megacam ugriyz + WIRCam Kg + 
Spitzer/IRAC 3.6/j,m4.5p,m photometry was used in the 
case of the two ECS candidates. 

Bouwens et al. (2015) derived a photometric redshift 
oi z = 7.00 for COSY-0237620370 over the CANDELS- 
COSMOS field and derived photometric redshifts of 8.1 
for the two sources over the CANDELS-EGS field (EGS- 
zs8-l and EGS-zs8-2, respectively), so the latter two can¬ 
didates were placed in the z ~ 8 sample of Bouwens et 
al. (2015). 

There was, however, some uncertainty as to both 
the robustness and also the precise redshifts of the 
GANDELS-EGS candidates from Bouwens et al. (2015). 
Prior to the present study, the use of the [3.6] — [4.5] color 
has never been systematically demonstrated to work for 
the identification of galaxies with redshifts of z > 7 de¬ 
spite there being ^5 prominent examples of z > 7 galax¬ 
ies with particularly red [3.6] — [4.5] colors (Bradley et 
al. 2008; Ono et al. 2012; Finkelstein et al. 2013; Tilvi 
et al. 2013; Laporte et al. 2014, 2015). Moreover, no 
Yios-band observations were available over either z > 7 
candidate from the GANDELS-EGS field in the Bouwens 
et al. (2015) selection to validate potential z > 7 galax¬ 


ies (though such observations have fortuitously become 
available as a result of observations made from the z9- 
GANDELS follow-up program [Bouwens et al. 2016]). 

The apparent magnitudes of the z = 7.1-8.5 galax¬ 
ies identified as part of the current selection are much 
brighter than the typical galaxy at z 8, as is evident in 
both the upper and lower panels in FigurelHl In fact, 3 of 
the sources from our current IRAG-red [3.6] — [4.5] > 0.5 
selection appear to represent the brightest z > 7.5 galax¬ 
ies known in the entire GANDELS program and consti¬ 
tute 3 of the 4 z ~ 8 candidates shown in the lower 
panel of Figure [H The only other especially bright 
Hieo,AB 25.0 z ^ 8 candidate shown in that lower 
panel is presented in the appendix (since it satisfies our 
[3.6] — [4.5] > 0.5 selection criteria using an independent 
set of photometry, i.e., Skelton et al. 2014). 

Interestingly enough, all 4 of the brightest candidates 
shown in the lower panel of Figure [9] are located in 
the GANDELS EGS field, providing a dramatic exam¬ 
ple of how substantial field-to-field variations in the sur¬ 
face densities of bright sources might be (though we note 
that EGSY-2008523660 is likely gravitationally lensed). 
This seems to be just a chance occurrence, as none of 
these candidates is clearly in a similar redshift window. 
The probability that the 4 brightest z ^ 8 sources in 
the GANDELS program would be found in the same 
GANDELS field (even if one is gravitationally lensed) 
is -1%[3 

There is only one source from our combined 2 ~ 8 selection 
with Bouwens et al. (2015) which would be potentially easier to 
select as a 2 : 8 galaxy over the CANDELS EGS field. It is 

presented in Appendix A. Its redshift is not well constrained (lying 
anywhere between 2 : 7.1 and 8.5), but it would be marginally 

easier to find over the CANDELS EGS field since the Bouwens et 
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Fig. 9.— {upper panel) Absolute magnitudes vs. redshift in the 
rest-frame UV for sources in our current photometric sample (solid 
red squares for the two sources from our sample with redshift mea¬ 
surements from spectroscopy and open red squares where the red¬ 
shift estimates derived from the photometry). For context, the 
absolute magnitudes and redshift measurements for other z > 6 
galaxies in the literature with spectroscopic redshift measurements 
from Lyo are shown (black squares: compiled from Jiang et al. 
2013, Finkelstein et al. 2013, Shibuya et al. 2012, Ono et al. 2012, 
and Vanzella et al. 2011). The gray-dashed line shows the evolu¬ 
tion of the characteristic magnitude M^y of the UV LF (Bouwens 
et al. 2015). Two sources from our sample, with redshift mea¬ 
surements from spectroscopy (the 4.7 (t one requires further con¬ 
firmation) are the brightest z > 7.5 galaxies discovered at such 
high redshifts and similarly for our bright photometric z ~ 8.6 
candidate (the downward arrow indicates the likely lensing mag¬ 
nification for this candidate: §3.4). {lower panel) Surface density 
of the full sample of z ~8 galaxies in the combined CANDELS 
and BoRG/HIPPIES fields (Bouwens et al. 2015, gray histogram). 
The shaded red squares indicate the position of our EGS-zs8-l and 
EGS-zs8-2 in the Bouwens et al. (2015) z ~ 8 selections, while 
the open green square indicates the position of the z ~ 8.6 candi¬ 
date EGSY-2008532660 (not identified as part of the Bouwens et 
al. 2015 z ~ 8 selection). Three of the sources from our selection 
represent the brightest-known galaxies at z > 7.5 (although one 
appears to be magnified from gravitational lensing). Interestingly 
enough, all three of the brightest, highest-redshift z ~8 candidates 
we have identified here (and four if one includes the source from 
Appendix A which is also in the Bouwens et al. 2015 z ~ 8 catalog) 
are located in only one of the CANDELS fields (CANDELS EGS), 
providing an example of how dramatic field-to-field variance can 
be for bright galaxies (see also Bouwens et al. 2015 and Bowler et 
al. 2015). 


al. (2015) z ~ 8 sample extends down to z ~ 7 over that field while 
the Bouwens et al. (2015) z ~ 8 samples over the other fields only 
extend down to z ~ 7.3. 


Previously, this point had been strongly made by 
Bouwens et al. (2015) in discussing the number of bright 
sources over the different CANDELS fields (Figure 14, 
Appendices E and F from Bouwens et al. 2015) and also 
quite strikingly by Bowler et al. (2015) in comparing the 
number of bright z ~ 6 galaxies over the UltraVISTA 
and UDS fields. 

6 . IMPLICATIONS FOR THE BRIGHT END OF THE Z ~ 8 

LF 

In this section, we will examine the implications of our 
present search results for the volume density of luminous 
galaxies in the z ~ 7-9 Universe. First, we estimate how 
complete we might expect our selection to be based on 
the [3.6] — [4.5] color distribution in fields where the red¬ 
shift can be constrained using deep U-band data (§ 6 . 1 ). 
Second, we make use of our search results and our com¬ 
pleteness estimates to set a constraint on the bright end 
of the z ~ 7-9 LF (§6.2). 

6 . 1 . [3.6]-[4--5] Color Distribution of z > 7 Galaxies and 
the Implications for the Completeness of our red 
IRAC Criteria and the [0111]+Hfl EWs 

In our attempts to identify bright z > 7 galaxies, we 
only consider those sources with red [3.6] — [4.5] > 0.5 
Spitzer/IRAC colors to ensure that the sources we select 
are robustly at z > 7 (see §3.2). However, by making this 
requirement, we potentially exclude those z > 7 galaxies 
which have bluer [3.6]-[4.5] colors, either due to lower- 
EW [OIII]-|-H/3 lines or simply as a result of noise in the 
photometry. 

To determine how important this effect is, we look at 
the [3.6]-[4.5] color distribution of galaxies which we can 
robustly place at a redshift z > 7 (where both lines in 
[OIII] doublet fall in the [4.5] band). The most relevant 
sources are those bright galaxies we can place at z > 7 
based on the available HSTj-ground-based photometry 
and which include deep flux measurements at l^m. Such 
measurements are available for the CANDELS GOODS- 
S, GOODS-N, UDS, and COSMOS fields, and a small 
fraction of the CANDELS EGS held. 

For our hducial results here, we only consider selected 
sources from those helds brightward of Hieo.AS = 26 
and with redshift estimates greater than z > 7.5. This 
is to ensure that we only include bona-hde z = 7.1-9.1 
galaxies (where the [OIII]-|-H/3 line in the 4.5^m band) in 
our selection. Photometric redshift errors often have an 
approximate size of Az~0.3 to this magnitude limit, and 
so to avoid z < 7 sources scattering into our selection, 
we kept our cuts fairly conservative. 

The list of such sources at such bright magnitudes is 
still somewhat limited at present, with only the bright 
^25.6-mag galaxy in the GANDELS South from Yan et 
al. (2012) and Oesch et al. (2012), a bright 25.7-mag 
galaxy in the GANDELS GOSMOS field from Bouwens 
et al. (2015), a bright ^25.5 mag source over the GAN¬ 
DELS GOODS-North held by Finkelstein et al. (2013), 
two bright sources over the GANDELS EGS held where 
Yi 05 -band photometry is available (EGS-zs 8 -l, EGS-zs 8 - 
2), and a third bright source over the GANDELS EGS 
held where the J 125 — Hieo color allows us to place it at 
z > 8 (EGSY-2008532660). 

Of these sources, hve out of the six sources have 
[3.6] —[4.5] colors in excess of 0.5, and therefore for sim- 
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TABLE 3 

Brightest z > 7.5 Galaxies over the CANDELS Fields and ~200 arcmin^ BoRG/HIPPIES area searched by 

Bouwens et al. (2015) 


ID 

R.A. 

Dec 

mAB‘" 

^phot 

[3.6] - [4.5] 

References* 

EGS-ZS8-1 

14:20:34.89 

53:00:15.35 

25.03±0.05 

7.9±0.4 

0.53T0.09 

[1], [8] 





(zspec=7.7302±0.0006) 


EGS-ZS8-2 

14:20:12.09 

53:00:26.97 

25.12±0.05 

7.6 ±0.3 

0.96T0.17 

[1] 





(2apec=7.4770±0.0008) 


EGSY-2008532660t 

14:20:08.50 

52:53:26.60 

25.26±0.09t 

8.57l“:l^t 

_Q ^QQ + 0.001n 

X'^spec —o.uoo_^ nna./ 

0.76T0.14 


GNDY-6379018085 

12:36:37.90 

62:18:08.50 

25.44±0.04 

7.508'= 

0 .88T0.11 

[6] 

BORGY-9469443552 

04:39:46.94 

-52:43:55.20 

25.56±0.20 

8.29l?:34 

— 

[1].[2].[3], [7] 

GSDY-2499348180 

03:32:49.93 

-27:48:18.00 

25.58±0.05 

7.84l°;i5 

0.08T0.09 

[1],[3],[4],[5] 

BORGY-6504943342 

14:36:50.49 

50:43:34.20 

25.69±0.08 

7 

< -"^^-3.17 

— 

[1] 

COSY-0235624462 

10:00:23.56 

02:24:46.20 

25.69±0.07 

'•o^_0.i8 

0.88T0.61 

[1] 

BORGY-2463351294 

22:02:46.33 

18:51:29.40 

25.78±0.15 

■7 Qq+0.59 

1 •»o_Q 21 

— 

[1] 

BORGY-2447150300 

10:32:44.71 

50:50:30.00 

25.91±0.20 

7 93+0.48 

1 •»0_Q jg 

— 

[1] 

BORGY-5550543040 

07:55:55.05 

30:43:04.00 

25.98±0.21 

' •^'^-5.63 

— 

[1] 

Median 





O.S2tlfo 



Other Possible Bright 2: > 7.5 Galaxies over CANDELS (Appendix A) 


EGSY-9597563148 

14:19:59.75 

52:56:31.40 

25.03±0.10 

8.19 

0.22±0.06 

[1] 


* References: [1] Bouwens et al. 2015, [2] Bradley et al. 2012, [3] McLure et al. 2013, [4] Yan et al. 2012, [5] Oesch et al. 
2012, [6] Finkelstein et al. 2013, [7] Schmidt et al. 2014, [8] Oesch et al. 2015b. 

^ The apparent magnitude of each source in the iLieo band. 

Maximum likelihood photometric redshift estimate from EAZY. 

^ Spectroscopic redshift determination (Finkelstein et al. 2013). Finkelstein et al. (2013) report a [3.6] — [4.5] color of 
0.98±0.14. 

t Photometric redshift estimate is not based on deep Y-band imaging observations over candidate of red [3.6] — [4.5] color 
as derived from the Bouwens et al. (2015) photometry. Inclusion in this list is based upon a red [3.6] — [4.5] color as 
derived from the Skelton et al. (2014) photometry. 

1 The flux of this source seems likely to be boosted by gravitational lensing (§3.4). 
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Fig. 10.— Range of Spitzer/IRAG [3.6] — [4.5] colors (blue his¬ 
togram) observed for bright z > 7 galaxies where we are confi¬ 
dent that the [OIII]-|-H/3 emission line falls in the [4.5] band (§6.1). 
Sources are included in this histogram if they are brighter than 
Hieo,AB = 26 and the redshift information available on these 
sources confidently place them at z > 7.5. This color distribu¬ 
tion is compared against the IRAC red [3.6] — [4.5] > 0.5 selection 
criteria we use (shaded red region). The sources presented here are 
the same as those sources in Table [3] with IRAC color measure¬ 
ments. Five out of 6 z > 7.5 galaxies show Spitzer/IRAG colors 
redder than 0.5. At face value, this suggests that our proposed 
[3.6] — [4.5] > 0.5 selection would identify 83(j(jg% of the bright 
2 > 7 population (but we emphasize this percentage is very uncer¬ 
tain due to the small numbers involved). The observed [3.6] — [4.5] 
color distribution also implies a minimum EW for [OIII]+H/3 of 
1300A. 


plicity, we will assume that our IRAC red selection is 
83% complete, but we emphasize that the completeness 
correction we derive from this selection is uncertain and 
could be much larger (as indeed one would expect if the 
[3.6] —[4.5] color measurement derived by Labbe et al. 
2013, i.e., ~0.4 mag, for the average stacked z ~ 8 galaxy 
are indicative). 

To investigate this possibility, we examined a slightly 
larger sample of objects over the four fields where we have 
photometric redshifts using K-band imaging. Consider¬ 
ing sources to a i?i 6 o-band magnitude limit of 26.2 over 
the CANDELS-UDS and COSMOS fields and 26.7 over 
the CANDELS GOODS-North and GOODS-South while 
extending the photometric redshift selection to z > 7.3, 
6 out of 9 sources satisfy the [3.6] —[4.5]>0.5 criterion. 
While this suggests the actual fraction of z > 7 galax¬ 
ies with such red IRAC colors may be less than 83%, 
this fainter sample is still consistent with our fiducial 
percentage. It also reassuring that our suggested selec¬ 
tion criteria would also apply to GN-108036, the bright 
JHi 4 o = 25.17 z = 7.213 galaxy found by Ono et al. 
(2012), given its measured [3.6] —[4.5] color of 0.58±0.18 
mag. 

We include a list of those sources and other sources in 
Table [3] from the Bouwens et al. (2015) catalog and also 
including the bright z = 7.508 galaxy from Finkelstein 
et al. (2013). In Figure [101 we present the [3.6]-[4.5] 
color distribution for the brightest sources we know to 
robustly lie at z > 7.5 based on spectroscopy or from 
the available HST-j-ground-based photometry for those 
sources that lie in regions of CANDELS with F-band 
observations or with J 125 — Tlieo colors red enough to 
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confidently place the sources at z > 8. 

The median [3.6] —[4.5] color that we measure is 
0.82 ^q'2o mag. Such a color implies a minimum EW 
of ^1300A for the [0III]+H/3 lines, assuming a flat stel¬ 
lar continuum and no line contribution to the [3.6] band. 
However, we emphasize that if there is also a substantial 
line contribution to the [3.6] band, e.g., from Hy, HJ, and 
[Oil], then the implied EW of the [OIII]+H/3 lines would 
be much larger. Eor example, adopting the line ratios 
from the Anders & Fritze-v. Alvensleben (2003), 0.2^0 
model would imply an EW of >2100A for the [OIII]-|-H/3 
lines. 

Some correction is required to the median [3.6] —[4.5] 
color measurement to account for the fact that the z > 7 
sources from the CANDELS EGS field were explicitly 
selected because of their red IRAC colors. If we assume 
that the intrinsic [3.6] —[4.5] color for sources over the 
CANDELS EGS held is ^0.6 mag (which is the value 
we hnd from the 3 candidates over the other CANDELS 
helds: see Table [3|) and the noise -I- scatter is ^0.4 (the 
value from the other helds), we compute a bias of 0.24 
mag from a simple Monte-Carlo simulation. Accounting 
for such biases reduces the median [3.6]-[4.5] color of the 
population by 0.24 mag, which implies a median EW of 
-800AandI500A, ignoring and accounting for a possible 
nebular contribution to the [3.6] band respectively. 

6.2. Volume Density of Bright z 8 Galaxies 

Here we use the search results from the previous section 
to set a constraint on the bright end of the z ^ 8 LF. 
We begin this section by calculating the total selection 
volume in which we would expect to hnd bright z > 7 
galaxies with our selection criteria. 

We will estimate the selection volumes in a similar way 
to the methodology used by Bouwens et al. (2015) in de¬ 
riving the LFs from the full CANDELS program. In 
short, we create mock catalogs over each search held, 
with sources distributed over a range in both redshift 
z ^ 6-10 and apparent magnitude (LIieo.AB = 24 to 
26). We then take the two-dimensional irys-band images 
of similar luminosity, randomly-selected z ^ 4 galaxies 
from the HUDF (Bouwens et al. 2007, 2011, 2015) and 
create mock images of the sources at higher-redshift us¬ 
ing the two dimensional pixel-by-pixel z ^ 4 galaxies as a 
guide (see Bouwens et al. 1998, 2003), adopting random 
orientations relative to their orientation in the HUDF 
and scaling their physical sizes as (I-|-z)“^’^, which is 
the approximate relationship that has been found com¬ 
paring the mean size of galaxies at hxed luminosity, as 
a function of redshift (Oesch et al. 2010; Grazian et al. 
2012; Ono et al. 2013; Holwerda et al. 2015; Kawamata 
et al. 2015; Shibuya et al. 2015). Individual sources were 
assigned UV colors based on their UV luminosity, using 
the 13-Muv relationship derived by Bouwens et al. (2014) 
and allowing for an intrinsic scatter ag of 0.35 at high lu¬ 
minosities {Muv.ab = [—22, —20]) as found by Bouwens 
et al. (2009, 2012) and systematically decreasing to 0.15 
at lower luminosities as found by Rogers et al. (2014). 

In addition to the HST images we created for individ¬ 
ual sources, we also constructed simulated ground-based 
and Spitzer/IRAC images for these sources which we 
added to the real ground-based -I- Spitzer/IRAC data. 
These simulated images were generated based on the 


mock 77i6o-band images we constructed for individual 
sources and convolving by the i7i6o-to-IRAC, Hieo-to- 
ground kernels that Mophongo (Labbe et al. 2010) de¬ 
rived from the observations. In producing simulated 
IRAC images for the mock sources, we assume a rest- 
frame EW of 300A for Ha-|-[NH] emission and 500A for 
[OHI]-|-H/3 emission over the entire range z = 4-9, a flat 
rest-frame optical color, and a TLieo-optical continuum 
color of 0.2-0.3 mag, to match the observational results 
of Shim et al. (2011), Stark et al. (2013), Gonzalez et 
al. (2012, 2014), Labbe et al. (2013), Smit et al. (2014, 
2015), and Oesch et al. (2013). 

We took the simulated images we created for individual 
sources and added them to the real HST, ground-based, 
and Spitzer/IRAC observations. These simulated images 
were, in turn, used to construct catalogs and our selec¬ 
tion criteria applied to the derived catalogs in exactly the 
same way as we applied these criteria to the real obser¬ 
vations (including excluding sources which violated our 
confusion criteria). 

Summing the results over all five CANDELS fields, 
we compute a total selection volume of 1.6x10® Mpc® 
per I-mag interval for galaxies with iLieo.AS magnitudes 
brightward of 25.5. If we assume that the present se¬ 
lection of z ^ 8 galaxies is complete, this would im¬ 
ply a volume density of <I.4xI0“® Mpc“® mag“^ and 
3.8l^;/xI0-® Mpc-® mag-^ for iLieo.AB 24.5-25.0 and 
Hieo^AB ^ 25.0-25.5 galaxies, respectively. We ignore 
the contribution of the z ^ 8.6 candidate galaxy EGSY- 
2008532660, given the evidence that it may be slightly 
magnified (§3.4). 

However, we cannot assume that the present selection 
of bright z ^ 8 galaxies is complete, since not every z ~ 8 
galaxy exhibits such a red [3.6] — [4.5] color. In the pre¬ 
vious section, we found that only 5 out of the 6 bright 
{Hieo.AB < 26), secure z >7 sources within CANDELS 
showed such red galaxies. Correcting the volume densi¬ 
ties given in the previous paragraph to account for this 
slight empirically-derived incompleteness (0.83to;i6), we 
estimate volume densities of <I.7xI0“® Mpc“® mag“^ 
and 4.7t2'7XlO“® Mpc“® mag“^ for ^ 24.5-25.0 

and Hieo,AB ^ 25.0-25.5 galaxies, respectively. The un¬ 
certainty in the completeness estimate is included in the 
error we quote for the volume density. 

It is interesting to compare these constraints on the 
volume density of bright z ^ 7.1-8.5 galaxies with other 
recent constraints which are available on the volume den¬ 
sity of bright end of the UV LFs at z ~ 4, z ^ 6, z ^ 8, 
and z ^ 10 galaxies from state-of-the-art studies (e.g., 
Bouwens et al. 2015). The results are shown in FigurefTIl 
and it is clear that our result lies somewhere midway be¬ 
tween the z ~ 7 and z ^ 8 LFs, as one might expect 
given the redshift distribution of the sources that make 
up our z = 7.1-8.5 sample. 

7. SUMMARY 

In this paper, we take advantage of the deep 
Spitzer/IRAC observations available over all five CAN- 
DELS in conjunction with the HST-|-ground-based data 
to conduct a search over 900 arcmin^ to find bright z ^ 8 
galaxies. To identify galaxies at such high redshifts, we 
select those galaxies with especially red Spitzer/IRAC 
[3.6] —[4.5] colors (i.e., >0.5), in the hopes of identifying 
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Fig. 11.— Estimated volume density of especially luminous 2 7.1-9.1 galaxies using the current IRAC-red search results over the entire 

CANDELS program (§6.2). This selection is assumed to be 17'^^% incomplete based on the results from §6.1 and excludes the 2 8.6 

candidate galaxy EGSY-2008532660 due to evidence that this galaxy may be magnified by a foreground source (§3.4). For context, the 
UV luminosity function results of Bouwens et al. (2015: B15) at 2 ~ 4, 2 ~ 6, 2 8, and 2 ~ 10 are also shown with the blue, cyan, red, 

and magenta open circles and solid lines, respectively. The present constraints on the UV LF are consistent with those previously derived 
by Bouwens et al. (2015). 


those z > 6 galaxies which show the presence of a strong 
[OIII]+H^ line in 4.5fim band. Such a selection is useful 
for the CANDELS program, given the lack of uniformly 
deep V-band observations over all five fields. 

Our selection yielded 4 z > 7 candidates brighter 
than an Hiqo,ab magnitude of 25.5. Each of these 
four selected candidates was requried to be undetected 
(<2.5cr) at optical wavelengths (<l^m), as defined by 
the inverse-variance-weighted mean flux measurement, 
be undetected in the Veoe-band (<1.5 (t), and show a 
hii — Ji 25 color redward of 1.5. 

Fortuitously, 3 of our 4 selected z > 7 candidates 
had deep F-band observations available from either 
deep ground-based observations or from the new z9- 
CANDELS follow-up program (Bouwens 2014; Bouwens 
et al. 2016) with HST. The available F-band observa¬ 
tions provide clear confirmation of the z > 7 redshifts we 
estimate for three of four candidates found in our search. 

The redshift estimates we obtain for three of our se¬ 
lected candidates lie significantly above z ~ 7, with 
EGS-zs8-2 having a redshift estimate of 7.6±0.3, EGS- 
zs 8 -l having a redshift estimate of 7.9±0.4, and EGSY- 
2008532660 having a redshift estimate of 8 . 6 to; 4 . 

We also obtained spectroscopic observations on two of 
our candidate z > 7 galaxies in the near-IR and find 
probable Lya lines in their spectra consistent with red- 
shifts of 7.4770±0.0008 and 7.7302±0.0006. The detec¬ 
tion of Lya emission for these candidates is significant 
at 4.7cr and 6.1 (t significance, respectively. The second 
of these sources was featured in Oesch et al. (2015b). 
Remarkably enough, a third candidate from our list was 
spectroscopically confirmed to lie at z = 8.683 by Zitrin 


et al. (2015). 

These sources represent the brightest z > 7.5 candi¬ 
dates we identified over the entire GANDELS program 
and are 0.5-mag brighter than z > 7.5 candidates identi¬ 
fied anywhere else on the sky. Goindentally enough, they 
all lie in the same GANDELS field, again suggesting large 
field-to-field variations for brightest z >7 galaxies. See 
also discussion in Bouwens et al. (2015) and Bowler et 
al. (2015). 

Using these candidates, we estimate the volume den¬ 
sity of bright {Hiqo,ab < 25.5) z >7 galaxies in the early 
Universe based on our selected sample and estimate that 
17^11% of z > 7 galaxies do not show such red colors. 
The volume density estimate we derive lies midway be¬ 
tween the volume density of luminous z ~ 6 galaxies 
Bouwens et al. (2015) derive and the volume density of 
luminous z ^ 8 galaxies. 

The median [3.6] —[4.5] color distribution for our selec¬ 
tion and other bright z > 7.5 galaxies from the literature 
is 0 . 82 to; 2 o (observed) and 0.581 q 2 o (correcting for 
the approximate selection bias: see §6.1). This strongly 
points to the existence of extremely high EW nebular 
emission lines in typical star-forming galaxies at z > 7. 
Assuming no contribution from nebular line emission to 
the [3.6] band implies a [OIII]-|-H/3 EW of ~800A. How¬ 
ever, allowing for contamination of the [3.6] band in ac¬ 
cordance with the expectations of Anders & Fritze-v. Al- 
vensleben (2003) would imply a median EW of ^1500 A. 
These results are in reasonable agreement but perhaps 
slightly higher than Smit et al. (2015) estimate for the 
IRAG-blue sources they selected at z = 6 . 6 - 6 .9. Smit et 
al. (2015) estimate a typical [OIII] -I- Hp EW of 1085A 
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for their selected sources. These estimates are similar 
albeit slightly higher than those estimated by Labbe et 
al. (2013), Laporte et al. (2014, 2015), and Huang et al. 
(2016). 

In the near future, we would expect the brightest z ^ 8 
galaxies to be identified within the ^1 deg^ wide-area 
UltraVISTA field (McCracken et al. 2012) by combin¬ 
ing the progressively deeper YJHKg observations with 
constraints from the optical Subaru-I-CFHT observations 
and Spitzer/IRAC observations from SPLASH (Capak et 
al. 2013) and SMUVS (Caputi et al. 2014). Another sig¬ 
nificant source of bright z 8 candidates will be the new 
BoRG[ 29 io] program (Trent! 2014), which uses a huge 
allotment of 500 orbits to cover a 500 arcmin^ area to 


> 26.5 mag depth (5 (t). 

We thank Robert Barone-Nugent, Daniel Schaerer and 
Dan Stark for valuable conversations. This work has 
benefited significantly from the public reductions of the 
SEDS program and hence the efforts of Matt Ashby, Gio¬ 
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also find a >3 (t line (putatively Lya) at 1.031^m (to ap¬ 
pear in D. Stark et al. 2016, in prep). We acknowl¬ 
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TABLE 4 

Additional bright (i^i6o,As < 25.5) z>7 source identified over the CANDELS fields 
UTILIZING THE SkELTON ET AL. (2014) CATALOGS FOR SOURCE SELECTION. 


ID 

R.A. 

Dec 

rriAB’^ 

[3.6] - [4.5] 

^phot^ 

References* 

EGSY-9597563148 

component-a 

component-b 

14:19:59.76 

14:19:59.78 

14:19:59.73 

52:56:31.40 

52:36:31.30 

52:36:31.70 

25.03±0.10 

25.73±0.14 

25.83±0.13 

0.53±0.26‘= 

S-lsto.s? 

[1] 


* References: [1] Bouwens et al. 2015 

^ The apparent magnitude of each source in the ffieo band. 

^ Maximum likelihood photometric redshift estimate from EAZY. 

^ This is the [3.6]-[4.5] color that Skelton et al. (2014) measure. Since our selection use the 
Spitzer/IRAC photometry from Bouwens et al. (2015) when available and Bouwens et al. (2015) mea¬ 
sure a [3.6]-[4.5] color of ~0.2 for the source, we do not include this source in our primary selection. 


^606 ^814 ^125 ^160 



Fig. 12. — HST/ACS 1606-^814 and HST/WFC3 Ji25Lfi60 postage stamp cut-outs {A" x 4^^) of one particularly bright (iLi60,AS < 25.5) 
z > 7 galaxy candidate selected by applying our IRAC-red criteria ([3.6] — [4.5] > 0.5) to the Skelton et al. (2014) photometric catalog. 
This source also has a red [3.6] — [4.5] color in the Bouwens et al. (2015) catalog, but do not quite satisfy our IRAC-red selection criterion. 
This galaxy appears to consist of two separate components (indicated with black hash marks and the labels “a” and “b” respectively). The 
photometry for each component is provided in Table [4| 


APPENDIX 

A. OTHER CANDIDATE Z >7 GALAXIES 

In addition to our application of our criteria to the catalogs Bouwens et al. (2015) compiled over a 750 arcmin^ 
search area within CANDELS, we also made use of the catalogs from the 3D-HST team (Skelton et al. 2014) over the 
same region. Our rationale to do so was to maximize the completeness of our selection for bright z > 7 galaxies. 

One additional z > 7 galaxy candidate is found which did not make it into our fiducial selection using the Bouwens 
et al. (2015) catalogs (because it had a measured [3.6] — [4.5] color of ~0.2 mag) IB We tabulate its coordinates, 
77i60,AB-band magnitude, [3.6] — [4.5] color, and estimated redshift in TablelH Postage images of the z ^ 8 candidate 
is presented in Figure [121 Model fits to the photometry Skelton et al. (2014) provide on the source, as well as the 
inferred redshift likelihood distribution, are also presented in Figure [131 

B. SOURCES USED TO VALIDATE OUR PROPOSED [3.6] - [4.5] > 0.5 SELECTION 

In §3.2, we considered a selection of sources from the four CANDELS fields with deep V-band observations to test 
the idea that we could use an IRAC color criterion, i.e., [3.6] — [4.5] > 0.5, combined with an optical dropout criterion 
to identify galaxies at z > 7 even in the absence of Y -band data. 

In Table m we provide a compilation of the 15 sources that we identified which satisfied the primary selection criteria 
fom the paper but which are brighter than 26.7 mag in the iLieo-band (and brighter than 26.5 mag over the CANDELS 
UDS and COSMOS fields). 


While such differences might seem to be a concern, the 
[3.6] —[4.5] colors we measure for the 4 other sources in our se¬ 


lection agree to <0.1 mag with the Skelton et al. (2014) values 
(~0.05 mag differences are typical). 
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Fig. 13. — Left: Best-fit SED models {blue line) to the observed HST -|- Spitzer/IRAC -h ground-based photometry {red points and error 
bars) for one especially bright (i^ieo AB < 25.5) z > 7 galaxies identified by applying our IRAC-red selection criteria ([3.6] — [4.5] > 0.5) 
to the Skelton et al. (2014) photometric samples. Right: The redshift likelihood distribution P{z) for the same candidate z >7 galaxy, as 
computed by EAZY using the Skelton et al. (2014) photometry. 


TABLE 5 

Sources in the 4 CANDELS fields with deep T-band data used to 

VALIDATE OUR PROPOSED SELECTION TECHNIQUE (§3.2) 


IDt 

R.A. 

Dec 

mieo^AB 

[3.6] - [4.5] 

^photX 

GNDY-6487514332 

12:36:48.752 

62:14:33.29 

26.4 

0 

'^•’^- 0.6 

7.66 

GNDY-7048017191 

12:37:04.805 

62:17:19.14 

26.2 

1 1 - 1 - 0 .3 
-*^•^- 0.2 

7.84 

GNWY-7379420231 

12:37:37.941 

62:20:23.14 

26.5 

n c-l-1-5 
^•^-0.5 

8.29 

GNWZ-7455218088 

12:37:45.529 

62:18:08.87 

26.5 

0 '7+0.2 
‘ - 0.2 

7.16 

GSDZ-2468850074 

03:32:46.889 

-27:50:07.45 

26.0 

1 1 +°-^ 
-*^•^- 0.1 

7.24 

GSWY-2249353259 

03:32:24.934 

-27:53:25.94 

26.1 

0 0 + 0-2 

8.11 

GSDY-2209651370 

03:32:20.964 

-27:51:37.02 

26.3 

^•'^- 0.8 

7.84 

COSY-0439027359 

10:00:43.90 

2:27:35.9 

26.6 

0 7 + 0.2 
‘ - 0.2 

7.33 

GOSZ-0237620370* 

10:00:23.76 

2:20:37.0 

25.1 

1 0 + 0-2 

7.14 

COSY-0235624462 

10:00:23.56 

2:24:46.2 

25.7 

'^•^- 0.1 

7.84 

UDSY-4133353345 

02:17:41.333 

-5:15:33.45 

25.8 

0 c+ 0-2 

U.O-o 2 

7.41 

UDSY-4308785165 

02:17:43.087 

-5:08:51.65 

26.3 

0 7 - 1 - 0.8 
'^•'-0.5 

7.84 

UDSZ-4199355469 

02:17:41.993 

-5:15:54.69 

26.5 

1 q + 2.6 

-*-•^- 0.6 

7.08 

UDSY-1765825082 

02:17:17.658 

-5:12:50.82 

26.3 

1 1 +°-^ 

7.93 

UDSY-5428621201 

02:16:54.286 

-5:12:12.01 

26.1 

-‘-•'^-0.7 

7.49 


Also featured in our primary selection. 

^ IDs are from the Bouwens et al. (2015) catalog. 

^ Photometric redshift computed based on the HST and deep ground-based 
data (i.e., not including constraints from the IRAC data). 


























